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Technical Memorandum X-53866 


Z TRANSFORM AND THE USE OF THE DIGITAL DIFFERENTIAL 
ANALYZER ASA PERIPHERAL DEVICE TO A 
GENERAL PURPOSE COMPUTER 


SUMMARY 


This work was undertaken to determine the feasibility of the Digital 
Differential Analyzer (DDA) as a peripheral device; that is, a computer system 
programmed in Fortran IV or some other standard language in which certain 
manipulations could be accomplished on the DDA by calling a subroutine. In 
this study, consideration was given to Large Scale Integrator Circuits (LSI) , 
because it was felt that suitable LSI would make the DDA a more feasible device. 

From the study, it appears that the DDA as a standard peripheral device 
with automatic patching is indeed feasible. However, implementation with LSI 
is questionable presently, but this may become feasible soon. 

A portion of time was allotted to the study of a control system analysis 
program using Z transform calculus. This work was performed to supplement 
the DDA work by extending the application of Z transform calculus to the 
analysis of the DDA. An analysis of IBM's control system analysis program 
was made, and several example problems were programmed using this program. 
The results and an explanation of these problems are given. 

The major advantage of the control system analysis program is that it 
can be used to derive difference equations, which are necessary to accomplish 
real-time flight simulation by digital computer techniques. It is felt that the 
program may be used in further analysis of the DDA. 



INTRODUCTION 


In the study of analog digital and hybrid computation, the author became 
interested in various approaches to the problem of solving differential equations. 
In reading about the Digital Differential Analyzer (DDA) , it occurred to the 
author that the speed of the DDA in making calculations would be a tremendous 
advantage if it could be used along with the digital computer, or possibly as a 
very good competitor with the general purpose digital computer. 

For the past year at Clemson University, the investigator has been 
teaching a course to seniors and first year graduate students in "Analog Digital 
and Hybrid Computation. " In this experience, many problems were encoun- 
tered; for example, patching of the Analog Computer, no hands on capability 
with digital, etc. From these problems, the feeling developed that the DDA 
may be suitable for alleviating some of the problems with computation. 

RESEARCH OBJECTIVES 


The following items were the objectives of this research study: 

1. To become familiar with the TRICE DDA. 

2. To search the literature to determine the state of the art of the DDA. 

3. To analyze the TRICE DDA using the Z transform calculus. 

4. To determine the feasibility of implementing a DDA using Large 
Scale Integrator (LSI) circuits. 

5. To define, design, and implement an interface for the DDA. 

6. To determine the numerical method best suited for the DDA and to 
develop the necessary software. 
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RESULTS OF STUDY 


Item 1 of the research objectives was accomplished partially by working 
on the hardware with Mr. Hugh Zeanah. A few problems were programmed 
and solved in order to obtain a basic working knowledge of the TRICE. More 
familiarity with the DDA was gained while studying the literature during the 
literature search, listed as item 2 of the research objectives. 

After obtaining a reasonable feeling for the state of the art, thought 
was given to some means of analyzing the TRICE DDA. IBM's control system 
analysis program ( 7090-MA-OIX) was already in existence and had previously 
been studied by Mr. Zeanah. This program has the capability to change 
control system relations from the S domain to the Z domain, and it was felt 
that this program would be helpful in analyzing the TRICE DDA, Therefore, 
effort was exerted to get the control system analysis program working to obtain 
the Z transform of several control problems. The results are given in the 
Appendix. 

Figure A-l of the Appendix shows the diagram of the system reduced to 
the S domain notation. Figure A-2 gives the necessary input data for obtaining 
a Z transform and root locus are given in Figure A -3. The printed results are 
presented in four different forms. The last three forms are variations of the 
format given in Reference 1, page 14. 

It is noticed that one value in the B matrix does not agree with that of 
the author's calculations. This was checked several times, however, and 
apparently there is a mistake in Reference 1. 

The first set of data is for 20 iterations, that is, the program carries 
out 41 terms (a=20 of the series for G* [ S] ) . 

G* (S) s -“r~ Yj g S+in-^- 
n=-a 

This relation represents the transfer function describing the output of a 
"sampler" that is sampling the impulse response of a continuous system G(S), 
where n is the number of samples and T is the time between samples. If a 
different number of terms is desired, the value of a is inserted in columns 
40 through 42 of the control card. The program automatically carries out 20 
iterations if columns 40 through 42 are left blank. 
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The second set of data is the same as the first except that the B matrix 
is left out. The B matrix consists of coefficients for the set of simultaneous 
equations used in calculating a Z transform. The B matrix is printed only if a 
-1 is set in field 17 (columns 49 through 51) of the control card. 

The third set of data is for 9 iterations (19 terms in the series for 
G* (S). The B matrix is also printed out. 

The fourth set of data is for 9 iterations. Here, only the Z transform, 
and not the root locus points from it, is computed. This is done by setting 
n=0 into field 16 (columns 46 through 48) where n should be equal to the degree 
of the resulting Z transform denominator polynomial. (For details of the 
program, see Reference 1. ) 

In order to further experiment with the control system analysis program, 
example problems were programmed [ 2, page 442] . The computer results 
compare favorably with Tou's results. (Refer to the Appendix. ) 

Harold Scofield of R-ASTR-FG was contacted concerning transforming 
control system problems. It was discovered that his group had experimented 
with obtaining Z transforms, but the group had not been able to carry the work 
as far as desired. It is believed that this program will be helpful in changing 
almost any of their control system problems into Z transform problems, and 
then into a difference equation to study the systems with a digital computer. It 
is suggested that more work be put into this effort, since the program will be 
helpful to R-ASTR-FG and will be a useful tool in analyzing the TRICE. 

In establishing the feasibility of implementing a DDA using LSI, it was 
necessary to get some idea of the state of the art of LSI. 

The present number of circuits obtainable from a pad leaves something 
to be desired. During fabrication, a portion of a pad may become defective and 
destroy circuit capability. If a portion of the circuit becomes defective there 
is presently no means of bypassing this defective portion, thus causing the 
entire board to become defective. It is also very difficult to build the LSI either 
automatically or by hand. Handwork is especially poor because of the small 
device involved. 

Generally, a custom integrated chip of LSI is very expensive but may 
become cheaper when the manufacturing techniques are improved. 

LSI is very suitable where a number of repeated circuits exist. 
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According to the Teledyne Systems Company of Northridge, California, 
a DDA implemented with LSI is feasible. They report of a Teledyne Electrically 
Alterable Digital Differential Analyzer (TEADDA), which is a parallel word, 
parallel computation Digital Differential Analyzer (DDA) with a speed of 
1 000 000 iterations (complete incremental solutions) per second and a reso- 
lution of one part per million. The TEADDA is programmed by interconnecting 
the computational elements electrically, controlled from a general purpose 
digital computer, tape reader, or other peripheral equipment. 

The TEADDA is implemented with Teledyne's hybrid LSI Mico-lElectronic 
Modular Assembly (MEMA) technique to achieve the high speed operation. High 
speed performance is achieved as a result of the small computer size and the 
MEMA being packaged so that short connections are used. The significant 
features of the TEADDA are as follows: 

1. 1-MHz iteration rate 

2. 20-bit resolution (one partner million) 

3. "Electronic Patchboard" for completely automatic programming 

4. 100 percent interconnectivity between 64 integrators (256 program- 
mable connections for each increment) 

5. Extrapolative trapezoidal integration 

6. Automatic register length selection to vary the incremental 
weighting over a range of 2 16 (approximately 65 000) 

7. Multiple incremental inputs to each integrator 

8. Ternary incremental communication 

9. Two's complement negative number representation 

10. MEMA packaging 

11. High reliability 

12. Low cost 

For more details on the TEADDA, see Reference 3, pages 161-169. 


5 



Items 5 and 6 of the objectives are being studied further. Also, more 
study is needed to firmly establish the feasibility of an LSI implemented DDA. 


CONCLUSIONS AND RECOMMENDATIONS 


It has been established that for some purposes, the DDA is more 
effective than a general purpose digital computer or an analog computer. It 
was discovered that digital computers and analog computers are extremely 
unwieldly in comparison with the DDA for some uses. In one case, it was found 
that solving an integral equation on a general purpose machine took 2 weeks, 
whereas on the DDA the time was 2 hours. 

There are real advantages to be gained in connecting the DDA to a high 
speed digital computer, since the integration operation on a DDA is comparable 
with the integration speed on a digital computer. 

When considering complexity and nonlinearity, the DDA is a very effective 
machine since the amount of equipment does not increase in proportion to the 
complexity of the problem. 

RECOMMENDATIONS FOR FURTHER STUDY 


It is recommended that the following further studies be performed to 
supplement and advance the present work. 

1. Continue to study control problems with the control system analysis 
program. With this program, it is expected that the Z transform of systems 
expressed in S transform can be derived. From the Z transform a difference 
equation can be derived that can be easily programmed on the digital computer. 
This technique can be used to further study the TRICE DDA. 

2. When sufficient information is obtained from item 1 above, modify 
the program or develop a program so that Z transforms can be automatically 
obtained from input and output pulse patterns. This would allow transforms to 
be derived for more complicated interconnections of TRICE elements. 

3. Make further studies of the availability and use of LSI to determine 
whether they are suitable for developing a new DDA. 

4. Continue work on items 5 and 6 of the Research Objectives. This 
should become a three or four man effort in order to be completed in about 1 year. 
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Figure A-l. System diagram, example 1. 

With the ratio of K /K held constant at 0. 6, the three runs will be 
r\ P 

made in this order: The first will be for the nominal values of Figure A-l, or 
perfect cancellation. The second will be for the case of the poles reduced in 
frequency by 20 percent. The third will be for the case with their frequency 
increased by 20 percent. Thus, for the first part of this triple run, with 
Kp = 1, having (arbitrarily) chosen to begin with the bottom feedback loop 

and work upward, we may work out an A -vector and M-vector sheet as 
follows! 1] : 
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Figure A -3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 1 of 7) 
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Figure A-3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 2 of 7) 
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Figure A-3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 3 of 7) 
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Figure A -3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 4 of 7) 
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Figure A -3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 5 of 7) 
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Figure A-3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 6 of 7) 
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Figure A-3. Root locus of a continuous system with 
calculated Z transform (20 iterations), (sheet 7 of 7) 



Shown in Figure A -4 is the block diagram of an error-sampled feedback 
control system. The sampling period is 0. 1 second. A zero order hold is used 
and the transfer function of the controlled system is 


G S ( S) 


30 

S (1 + 0. IS) 



Figure A -4. System description, example 2. 
The forward transfer function is 


G(S) 


30 (1 
S z (1 + 0. IS) 


-TS 

e ) 


The Z transform associated with G(S) is found to be 


G(Z) 


1. 104 (Z + 0.718) 
(Z-l) (Z - 0.368) 


The computer results (Fig. A -5) compare identically. 


G(Z) 


1. 1036365 Z + 0. 7972486 
Z 2 -l. 3678794 Z + 0.36787944 
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Card # 


Input Data Card 


1 

m 

10 

1 

0 

-11 

3 

2 

1 

D 

■ 

1 

1 

1 

■ 

■ 

-2 

■ 


2 


. 5 

. 5 

2 

• 

0. 

6. 

12. 

. i 



3 

10. 

m 

■ 

■ 

1. 

300. 

B 

B 

0. 

0. 

4 

10 


i. 

0. 

1. 






5 

2 

1 

-i 

0 

2 

0 

-1 

2 

3 

-1 

■ 

■ 

■ 

■ 

■ 

■ 

■ 



Computer Results For Example 2 


FACTORED DENOMINATOR OF Z TRANSFORM 


-. 36787944E 

00 

. 10000000E 

01 

-. 10000000E 

01 

. 10000000E 

01 

Z TRANSFORM NUMERATOR 




. 79272486E 

00 

. 1 1036365E 

01 

. 00000000E 

00 




POLYNOMIAL DENOMINATOR OF Z TRANSFORM 

. 36787944E 00 13678794E 01 

. 10000000E 01 

Figure A -5. Data for example 2. 

The following example demonstrates the ability of the control system 
analysis program to transform a ratio of polynomial, part of which may be in 
S and part in Z, into the Z domain. The block diagram of the system is given 
in Figure A -6. Note that this is Example 2 modified by additional networks 
that are already expressed in Z notation. 
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( Z - 0. 5) (Z- 0. 333) 


-TS 

r 1 " e 


30 

Z ( Z + 0. 5) 


G h ^ S 


S(1 + 0. 15) 


Figure A-6. System description, example 3. 

30 (l-e~ TS ) (Z - 0.5) (Z - 0.333) 

S(S) (1 + 0. 15) Z(Z + 0. 5) 

300 (l-Z" 1 ) (Z - 0. 5) (Z - 0. 333) 

S(S) (S + 10)(Z) (Z + 0. 5) 


300 (Z-l) / Z (Z - 0.5) (Z - 0.333) 
Z (Z + 0. 5) (S*) (S + 10) 


The zero order hold is not considered as part of the system poles. 
Hence the system poles are: 

(Z+0) (Z+0.5)(Z-1) (S+10) 

1 


The pole at the origin, S, is entered as Z-l. 

The parameters are now entered for the entire system expression: 

G = (Z-l) ( Z-0. 5) (Z -0.333) (300) 

(Z a + 0. 5Z z + 0Z+0) (S 3 +10S 2 +0S+0) 


From the two equations immediately above, the input data are obtained 
and are given in Figure A-7. 
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14 
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Computer Results For Example 3 


Factored Denominator of Z Transform 


36787944E00 
10000000E01 
. 50000000E00 
. 00000000E00 


. 10000000E01 
. 10000000E01 
. 10000000E01 
. 10000000E01 


Z Transform Numerator 


. 13199240E00 47659026E00 

12659967E00 . 11036323E01 

. 00000000E00 


Polynomial Denominator of Z Transform 
. OOOOOOOOEOO . 18393972E00 

31606028E00 86787944E00 

. 10000000E01 


Figure A -7. Data for example 3. 
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The computer results may be checked by considering the expression in 
the form 

( _TS\ 

1-e l 

U ~ Z (Z+0.5) (S 2 ) (S+10) ’ 


which transforms into 


G (Z ) 


(Z-0. 5) (Z-0. 333) (1.104) (Z + -0.71 8) 
Z (Z+0.5) (Z-l) (Z-0. 368) 


and when multiplied out becomes: 


G(Z ) 


1. 104 Z 3 - 0. 127 Z 2 - 0. 476 Z + 0. 132 
Z 4 - 0. 868 Z 3 - 6. HeKTz 2 + 0. 1840 Z 


This checks with the computer results given in Figure A -7. 

When rounding off the computer results it is observed that they compare 
with the calculated values. 

0, 132 - 0. 477 Z - 0. 127 Z 2 +1. 104 Z 3 = , 

0. 184 Z-0. 316 Z 2 - 0. 868 Z 3 + Z 4 


20 



REFERENCES 


1. IBM Application Program. H20-0052-1 Program Reference Manual 
7090/7094 Control System Analysis (7090-MA-01X) . 

2. Tou, Julius T. : Digital and Sampled Data Control Systems. McGraw-Hill 
Book Co. , Inc. , N. Y. , Toronto, London, 1959. 

3. Hyatt, G. P. and Ohlberg, G. : Electrically Alterable Digital Differential 
Analyzer. Teledyne Systems Co. Northridge, Calif. , Spring Joint 
Computer Conference, American Federation of Information Processing 
Societies, N. Y. , N. Y. , 1968, p.*l61. 


21 



BIBLIOGRAPHY 


Burtnelt, Thomas W. : A Microcircuit Integrator for a Digital Differential 
Analyzer. Report No. 1-65-36. 

Daniel, James W. : A Partial Error Analysis for the Solution of Differential 
Equations in Simulation: A look at Fowler's Z Transform Root-Locus Method. 
IBM Journal, Nov. 1966, pp. 472-476. 

Eckes, Harry R. : Disign, Test, and Application of a High-Speed Iterative 
Analyzer. University of Arizona N67-31232, Tucson, Arizona, NASA. 

Extension of the Application of the Z Transform Calculus to Analysis of 
Incremental Computers (TRICE DDA) . Accession No. N68-10205, NASA CR or 
TMX or AD No. CR-65-787. 

Fowler, M. E. : A New Numerical Method for Simulation. IBM Systems 
Research and Development Center, Palo Alto, Calif. , Simulation Councils, Inc. , 
Vol. 4, No. 5, May 1965. 

Fowler, M. E. andWarten, R. M. : A Numerical Integration Technique for 
Ordinary Differential Equations with Widely Separated Eigenvalues. IBM 
Journal, Sept. 1967, pp. 537-543. 

Goldman, M. W. : Design of a High Speed DDA. Proceedings - Fall Joint 
Computer Conference, Martin Company, Baltimore, Md. , 1965. 

LeWallen, R. A. and Trauboth, H. H. : The Use of DSL/90 for the Simulation 
of Stabilized Platform Dynamics. NASA TMX-53628, Marshall Space Flight 
Center, Huntsville, Alabama. 

Mayorov, F. V. : Electronic Digital Integrating Computers, Digital Differential 
Analyzers. London Iliffe Books Ltd. , New York American Flsevier Publishing 
Co. , Inc. Originally published in the USSR under the title Elektronnye 
Tsifrorye Integriruynshchie Mashiny, Translated by Scripta Technica, Inc. 

Mitchell, J. R. ; Moor, J. W. ; and Trauboth, H. H. : Digital Simulation of an 
Aerospace Vehicle. Proceedings - ACM National Meeting, Marshall Space 
Flight Center, Huntsville, Alabama, 1967. 


22 



BIBLIOGRAPHY (Concluded) 


Numerical Techniques for Real-time Digital Flight Simulation. IBM Manual 
E20-0029, 

O'Grady, E. P. : A Hybrid-Code Differential Analyzer. N67-13605, NASA, 

May 1965. 

Owen, P. L. ; Partridge, M. P. ; Sizer, M. A. ; and Sizer, T. R. H. : Corsair, 

A Digital Differential Analyzer. TN No. I. A. P 1123, December, 1960. 

Scheiman, G. : DDA Organization for an Analytic Platform Transformation 
Computer. Army Missile Command, Huntsville, Ala. , January 16, 1967. 

Scheiman, Gerald: Digital Differential Analyzer Organization for an Analytic 
Platform Transformation Computer. Report No. R. G. TR67-2, U. S. Army 
Missile Command, Redstone Arsenal. Ala. 

Smith, M. G. and Notz, W. A. : Large Scale Integration from the Users Point 
of View. Fall Joint Computer Conference, 1967, pp. 87-94. 

Starr, David A. and Johnson. Jens J. : The Design of an Automatic Patching 
System, Simulation Councils, Inc,, Vol. 10, No. 6, June 1968. 

Technique for Solving Nonlinear Differential Equations. A FIPS Conference 
Proceedings, Vol. 25, Spring Joint Computer Conference, 1964, pp. 169-179. 

Trauboth, H. H. : Digital Simulation of General Control Systems. Simulation 
Councils, Inc. , June 1967. 

TRICE Lecture Notes CSP156. Packard Bell Computer, Los Angeles, Calif. 

Wait, John V. : A Hybrid Analog-Digital Differential Analyzer System. Uni- 
versity of Arizona N66-1 1389, Tucson, Arizona, NASA. 

Wodke, Hank: TRICE Programming and Operations Course Notes. Raytheon Co. 

Zimmerle, D. F. : Development of an Improved Integration Scheme for Digital 
Differential Analyzers. NCR Report 67-3-10M, Contract NAS8-21213, Monthly 
Technical Progress Report to Marshall Space Flight Center, NASA, Huntsville, 
Ala. , May 1, 1968. 


23 



APPROVAL 


TMX-53866 


Z TRANSFORM AND THE USE OF THE DIGITAL DIFFERENTIAL 
ANALYZER AS A PER I PHERAL DEV I CE TO A 
GENERAL PURPOSE COMPUTER 


By Donald Nalley 


The information in this report has been reviewed for security classifi- 
cation. Review of any information concerning Department of Defense or Atomic 
Energy Commission programs has been made by the MSFC Security Classifica- 
tion Officer. This report, in its entirety, has been determined to be unclassi- 
fied. 

This document has also been reviewed and approved for technical 
accuracy. 



Hugh w.^zeanah 

Section Chief, Digital 


raphics and Simulation Section 



Director, Computation Laboratory 


24 


MSFC— R? A, Ala 



DISTRIBUTION 


TMX-53866 


Internal 

DIR 

Dr. von Braun 
AD-S 

Dr. Stuhlinger 

DEP-T 

Dr. E. Rees 

S&E-DIR 

Mr. H. Weidner 

PD -DIR 

Dr. W. R. Lucas 
Dr. W. A. Mrazek 
Mr. F. L. Williams 

PM-PR-M 

S&E-COMP-DIR 
Dr. H. Hoelzer 
Mr. Carl Prince 

S&E-COMP-C 
Dr. H. Kerner 

S&E-C OMP-CS 

Dr. H. Trauboth (3) 
Mr. D. Marion 

S&E-C OMP-CR 
Mr. A. Dean 


S&E-CSE-DIR 

Dr. W. Heaussermann 

S&E-ASTR-DIR 
Mr. B. Moore 

S&E-ASTN-DIR 

Mr. K. L. Heim burg 
Dr. R. Head 

S&E-ME-DIR 
Dr. M. Siebel 

S&E-SSL-DIR 

Mr. G. B. Heller 

S&E-QUAL-DIR 
Mr. Grau 

S&E-R-DIR 

Mr. W. G. Johnson 

A&TS-MS-IL (8) 

A&TS-MS-IP (2) 

A&TS-MS-H 

A&TS-TU (6) 

A&TS-PAT 

Mr. L. D. Wofford, Jr. 
External 


S&E-COMP-S 
Dr. Polstorff 
Mr. R. Lawrence 

S&E-AERO-DIR 
Dr. E. Geissler 


Dr. R. Barfield 

Dept, of Mechanical Systems 

Engineering 

Univ. of Alabama 

University Ala. 35486 


25 



DISTRIBUTION (Concluded) 


TMX-53866 


External (Concluded) 

D. W. N alley (15) 

Dept, of Electrical Engineering 
Clemson University 
Clemson, S. C. 29631 

Scientific and Technical 
Information Facility (25) 

P. O. Box 33 

College Park, Maryland 20740 

Attn: NASA Representative, S-AK/RKT 


26 



